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Nanoscale electronic inhomogeneity in optimally doped single-layer Bi2Sri.6Lo.4Cu06+a, {L- 
Bi2201, L—Lsl and Gd) as well as bilayer Bi2Sr2CaCu208+« (Bi2212), superconductors has been 
investigated by the low-temperature scanning-tunneling-microscopy and spectroscopy. The inhomo- 
geneous gap structures are observed above and below Tc on Gd-Bi2201. Increasing the out-of-plane 
disorder which is controlled by mismatch of the ionic radius between Sr and L results in the enhance- 
ment of pseudogap region with suppressed coherence peaks in the conductance spectrum, making 
the gap average larger and the gap distribution wider. The degree of electronic inhomogeneity 
measured by the distribution width of the gap amplitude is shown to have a correlation with Tc in 
optimally doped Bi2201 and Bi2212. 

PACS numbers: 74.25.Jb, 68.37.Ef, 74.72.Hs 



To explore the mechanism of high-Tc superconduc- 
tors (HTSC's), tunneling spectroscopy is a powerful tool 
since it provides direct information on the electronic den- 
sity of state. Especially, recent experimental results of 
scanning tunneling microscopy and spectroscopy (STM 
and STS) have revealed the nanoscale electronic inho- 
mogeneity which shows up as spatial modulation of gap 
amplitudeiiiSiSi^ The coherence peaks on both side of 
the gap edge are suppressed as the gap amplitude gets 
larger, which is similar to the pseudogap in normal state 
of underdoped cuprates^ The inhomogeneous electronic 
structure is suggestive of quite unconventional supercon- 
ducting state and gives a strong impact on the explo- 
ration of the mechanism of high-Tc superconductivity. 
Disorder introduced by the chemical doping processes is 
a possible origin of the electronic inhomogeneity^ In fact, 
for Bi2Sr2CaCu208+5 (Bi2212), the doping is controlled 
by the amount of excess oxygen S near the BiO layers 
whch is most likely a source of the microscopic electronic 
inhomogeneity 

In addition to the oxygen nonstoichometry, the disor- 
der introduced by cationic substitution is also expected 
to play an important role in HTSC. It has been suggested 
that disorder outside the Cu02 planes might be responsi- 
ble for the large difference in Tc, ranging from 30 to 90 K, 
among different cuprates with a single Cu02 plane in the 
unit cell. Recently, Fujita et al. have demonstrated that 
the cationic disorder in the layer next to the Cu02 plane 
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in single-layer cuprates substantially reduces Tc without 
changing the carrier density^ However, at present it is 
not clear how cationic disorder outside the Cu02 planes 
(out-of-plane disorder) is related to the microscopic elec- 
tronic inhomogeneity. 

In this paper, we report on a direct proof of the system- 
atic variation of nanoscale electronic inhomogeneity ob- 
served by the STM and STS measurements on rare-earth- 
substituted Bi2Sr2Cu06-i-5 (L-Bi2201) with the same 
doping level, which were well characterized in the preced- 
ing worki. The La or Gd ions are partially replaced for 
the Sr ion in the SrO plane, next to the Cu02 plane. The 
ionic radius of L (1.22 A for La.^+ and 1.10 A for Gd^+) 
is smaller than that of Sr^+ ion (1.27 A), so that this mis- 
match introduces structural disorder in the SrO planes. 
Since the mismatch (and hence disorder) is larger in Gd- 
Bi2201, a comparison between La and Gd would sort out 
the correlation between out-of-plane disorder and elec- 
tronic inhomogeneity in the Cu02 planei^ We show that 
the nanoscale gap inhomogeneity is strongly correlated 
with this out-of-plane disorder and the degree of inho- 
mogeneity evaluated by gap amplitude distribution is in- 
timately related to the superconducting transition tem- 
perature Tc- 

The L-Bi2201 (nominally Bi2Sri.6Lo.4Cu06+5) single 
crystals were fabricated by the traveling-solvent-floating- 
zone (TSFZ) methodpL Tc is 34 K for La-Bi2201 and 14 
K for Gd-Bi2201. The samples were prepared to have op- 
timal doping level by oxygen annealing under the same 
condition (at 650 °C and for 48 h). The constant carrier 
density is confirmed by the measurement of the trans- 
port properties, the slope of T-linear resistivity, and the 
magnitude of Hall coeflicientii For comparison, the op- 
timally doped Bi2212 (nominally Bi2.iSri.gCaCu208+5, 
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Tc=93 K) was also examined. A homebuilt STM unit in- 
stalled in a low-temperature chamber cooled by liquid He 
has been used for the present measurements. The tem- 
perature of the sample stage was typically fixed at 9 K 
except for temperature dependent measurements on Gd- 
Bi2201. The single crystals were transferred to the STM 
unit soon after the cleavage inside the low-temperature 
chamber under the ultrahigh-vacuum atmosphere. The 
Pt/Ir tip was cleaned by a high- voltage field emission 
process with Au film target just prior to the STM and 
STS observations. 

Figure 1 shows STM-STS mcsurements of temperature 
dependence on Gd-Bi2201 below and above Tc (at 5 K, 
9 K, and 26 K). In Figs. l(a)-l(c), we show (a) to- 
pography, (b) the spatial distribution of gap (A map) 
and (c) the spatial variation of the conductance spectra 
[G{x, V) = dl/dV] along the white line in Fig. 1(b) mea- 
sured on Gd-Bi2201 at 5 K (< Tc). The gap A was de- 
fined as half of the maximum-conductance peak-to-peak 
voltage in each spectrum. In some spectra, the gap peak 
is suppressed in the occupied-state side (negative-bias 
side). In these cases, A is defined by the peak voltage 
in the unoccupied-state side. Figures l(d)-l(f) show (d) 
topography, (e) A map, and (f) spatial variation of the 
conductance spectra at 9 K (< Tc). At both 5 K and 9 
K, the spatial inhomogeneous gap structures were clearly 
observed basically similar to these of Bi2212 as observed 
in previous workaiiSii^ and our experimental results [see 
Fig. 2(d) and 2(g)]. The comparison between the differ- 
ent samples will be discussed later. At this stage, it is im- 
portant to check the influence of the thermal broadening 
effect on the conductance spectra based on the tempera- 
ture dependence measurement, since the temperature of 
material dependent measurements discussed later (9 K) 
is close to Tc (=14 K) in the case of Gd-Bi2201. However, 
as easily recognized by Figs. 1(c) and 1(f), the broadened 
gap peaks and the distribution of A (from 15 meV to 120 
mcV) at 9 K are quite similar to those of 5 K. Based on 
this fact, we conclude that the broadened spectra of Gd- 
Bi2201 were not due to the thermal effect but to some 
intrinsic characteristics of the electronic states such as 
the lifetime broadening due to the scattering by the dis- 
order. Moreover, we should note that the inhomogeneous 
gap structures survive even above Tc. In Figs. l(g)-l(i), 
we show (g) topography, (h) A map, and (i) spatial vari- 
ation of the conductance spectra at 26 K (> Tc). The 
spectral features above Tc and the spatial extent of the 
inhomogeneity of the A are again similar to those ob- 
served at 5 K and 9 K. These results indicate that the 
inhomogeneous electronic structure and the conductance 
spectra are relatively insensitive to the temperature up 
to 26 K even to the onset of superconductivity. This im- 
plies that most of the gap structures observed here can 
be attributed to the pseudogap origin rather than the su- 
perconducting origin. The existence of the gap structure 
above Tc is consistent with the previous results of STS^S 
angle-resoleved photoemission spectroscopjiiS (ARPES) 
and NMR (Ref.ii) experiments on Bi2201. 
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FIG. 1: (Color online) The temperature dependence of STM and 
STS measurement on Gd-Bi2201. (a) Topography at 5 K {It = 0.05 
nA, Vtias = 0.4 V), (b) A map at 5 K on the same area of (a) 
(A= 52 meV, (ta > 22 meV within the field of view), (c) The 
raw conductance spectra G{x,V) along the white line in (b). (d) 
Topography at 9 K {It = 0.1 nA, Vi,ias = 0.1 V). (e) A map at 9 K, 

(f) The raw conductance spectra G{x, V) along the white line in (e). 

(g) Topography at 26 K {It = 0.1 nA, Vuaa = 0.1 V). (h) A map 
at 26 K. (i) The raw conductance spectra G{x, V) along the white 
line in (h). The length-scale bars are common with topographies 
and zA maps. The color scales of A are common with three A 
maps. The magnitude of conductance is normalized by G{x, O.IV). 



From our experimental results of the A inhomogeneity 
above Tc, it is necessary to reconsider the meaning of the 
A inhomogeneity including that below Tc. The degree of 
A inhomogeneity (for example, the standard deviation of 
A) may not be directly attributed to superconductivity. 
However, the A inhomogeneity generated by pseudogap 
may play an important role as the background of super- 
conductivity and has some correlation with Tc. Anyway, 
at present, although the meaning of the degree of A inho- 
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mogcncity is ambiguous, it has some significant features 
for high Tc superconductivity. The relation between the 
deviation of A and measured by microscopic probe 
of STM and STS is considered to be corresponding to 
the relation between the mismatch of ion radius and 
measured by macroscopic method^ 




FIG. 2: (Color online) (a)-(c) Topographic images on (a) Bi-2212 
{It = 0.1 nA, Vtias = 0.07 V), (b) La-Bi2201 (h = 0.1 nA, Vbias = 
0.1 V), and (c) Gd-Bi2201 (h =0.1 nA, Vbias = 0.1 V) at 9 
K. (d)-(f) A maps of 20 nm X 20 nm area on (d) Bi2212. (e) La- 
Bi2201. (f) Gd-Bi2201 at 9 K, including the topography area in (a), 
(b), and (c), respectively. The length-scale bars are common with 
topographies and A maps, (g)-(i) The raw conductance spectra 
G{x, V) along the white line in each A maps at 9 K. The magnitude 
of conductance is normalized by G(x,0.1V). 



Figures 2(a)-2(c) show the topographic images mea- 
sured on (a) Bi2212, (b) La-Bi2201, and (c) Gd-Bi2201 



at 9 K, respectively. In each topography, supermodula- 
tion in the BiO layer is observed. The whole area in each 
sample is covered with these supermodulation structures 
with periodic length of 2.5±0.4 nm at least within an 
area of 70 nm x 70 nm, confirmed by the observation 
of wide-area topography. The condition of the exposed 
BiO layer is considered to have no drastic difference be- 
tween Bi2212 and L-Bi2201. Figures 2 (d)-(f) are the 
A maps measured on (d) Bi2212, (e) La-Bi2201 and (f) 
Gd-Bi2201 within the area of 20 nm x 20 nm. Figures 
2(g)- 2(i) are spatial variations of the conductance spec- 
tra measured along the white line in each A map. In the 
gap map of Bi2212, one can easily see that A varies on a 
nanometer scale as observed in previous works. The spa- 
tial distributions of A on La- and Gd-Bi2201 are basically 
similar to that of Bi2212 but the distribution of the gap 
amplitude is much broader. It is seen that, as disorder 
becomes stronger, in the order of Bi2212, La-Bi2201, and 
Gd-Bi2201, the fraction of large-gap area (dark colored 
area) increases. The length scale of the inhomogeneity 
(2^a) defined as the full width of half maximum of au- 
tocorrelation analysis is 2^a —2.7 nm for La-Bi2201, 2.2 
nm for Gd-Bi2201, and 3.1 nm for Bi2212. It is slightly 
smaller in Gd-Bi2201. These lengths arc comparable with 
the superconducting coherence length of the order of 
a few nm. Note that, if follows the BCS relationship 
^0 = Hvp/ttAo (X l/Tc, of Gd-Bi2201 should be longer 
by a factor of 2 or more than that of La-Bi2201. There- 
fore, ^A does not scale with Tc, but it is roughly the same 
among Bi-based cupratcs. 

To see the evolution of the spectrum with increasing 
disorder the conductance spectra are shown in Figs. 3(a)- 
3(c). Each curve is obtained by averaging all the spectra 
with identical gap amplitude within the ± 2 meV in each 
sample. Compared with the spectra of Bi2212, the spec- 
tra of L-Bi2201 exhibit following distinct features. 

(1) A is more widely distributed in the range between 
10 and 100 meV, and the larger gap region (A >50 meV) 
dominates. For large gaps, the spectrum has V-shaped 
structure without showing sharp coherence peaks. 

(2) The zero-bias conductance (ZBC) is appreciably 
large and increases as A decreases. The ZBC in Gd- 
Bi2201 is larger than that in La-Bi2201. 

The distributions of A at 9 K are shown in histograms 
in Fig. 3(d), Bi2212, (e), La-Bi2201, and (f), Gd-Bi2201. 
One can see that the distributions of A in La- and Gd- 
Bi2201 are much broader, and show a long tail in the 
high-energy side. The parameters that characterize the 
A distribution (A and a a, the spatial averaged gap and 
the standard deviation, respectively) arc A— 42 meV and 
(TA=7 meV for Bi_2212, A=36 meV and crA = 17 mcV for 
La-Bi2201, and A=48 meV and crA=22 meV for Gd- 
Bi2201. If the inhomogeneity of A is simply generated 
by the spatial variation of local carrier density, the av- 
erage gap amplitude A should be identical between La- 
and Gd-substituted systems, since they have the same 
doping levelii However, in Gd-Bi2201, A is larger than 
that in La-Bi2201, reflecting the fact that the fraction of 



4 



3 ■ (a) 




-151) -100 -50 — — So — anso « 2" « " ™ 1«» ™ 



I 







L3-BI2201 




4. 




■ — ia.iiiiiV 






' — ^ 









Gd-Bi2201 SK 




\ 1 i' /// 



150 -100 -50 50 100 150 
V....,,.(n,V) 




FIG. 3: (Color online) (a)-(c) The A-resolved spectra on (a) 
Bi2212, (b) La-Bi2201, and (c) Gd-Bi2201 at 9 K. The spectra 
were obtained by collecting and averaging from all the raw spectra 
data within the ± 2 mV of each A, except for the case of A > 80 
meV within the it 4 meV. The black dot in each zero-bias posi- 
tion indicates the average of zero-bias conductance, (d)-(f) The A 
histogram on (d) Bi2212, (e) La-Bi2201, and (f) Gd-Bi2201 at 9 K. 



large gap region increases in the former. This indicates 
that disorder in the SrO layers influences the A distri- 
bution and the average value A in Bi2201 and that A 
is not directly related to Tc nor to the total carrier den- 
sity. However, it is still an open question whether or not 
disorder gives rise to modulation of local carrier density. 

Note that the conductance spectra with identical gap 
amplitude are significantly different between La-Bi2201 
and Gd-Bi2201. This is in contrast to the case of 
Bi2212 with different doping levels;^ where the conduc- 
tance spectrum at the position of identical A is found 
to be nearly the same for all doping levels. This differ- 
ence as well as the appreciable ZBC may arise from much 
stronger disorder in Bi2201. The large-gap spectrum 
with V-shaped structure is similar to that observed for 
heavily underdoped Bi2212(Refs>ii^) and lightly doped 
Naa;Ca2-2:Cu02Cl2 (RefsiiSii^), although the spectrum 
shape is slightly different from that for underdoped 
cuprates. possibly due to different location of the Van 
Hove singularity. The different spectral features includ- 
ing ZBC in Bi2201 will be discussed in more details 
elsewhere 

We see that A of Bi2212 (A=:42 meV) is higher than 
that of La-Bi_2201(A=:36 meV) while lower than that of 
Gd-Bi2201 (A=48 meV). There is no apparent correla- 



tion between A and Tc- On the other hand, Tc is found 
to be scale with a a, as displayed in Fig. 4, in which Tc is 
plotted as a function of a a including the available data 
for Bi2212 near optimal doping,ii^ii^ For Bi2212, we 
plot a few data measured on nominally overdoped sam- 
ples, but their Tc values are so high (Tc > 82 K) that their 
doping level is considered to be near optimal. Tc is found 
to decrease linearly with of cta over a wide range of Tc 
values, irrespective of the number of the number of Cu02 
planes and of the dominant disorder sites. In Bi2212, dis- 
order is presumably excess oxygen atoms near the BiO 
layers, but cationic disorder also exists in the SrO lay- 
ers in the form of Bi2-|-a:Sr2— a; CaCu208+5. The resent 
STM-STS observation by McElroy et alM. suggests that 
the excess oxygen atoms are a dominant source of gap 
inhomogeneity in Bi2212. On the other hand, cationic 
disorder in SrO layers dominates in Bi2201 and is prob- 
ably responsible for much broader gap distribution and 
maybe for lower Tc 

The inset of Fig. 4 illustrates a schematic variation 
of Tc as functions of doping and disorder (the present 
experimental conditions are indicated by open circles). 
We have shown that the spatial distribution of A as well 
as the averaged A is largely different among the sam- 
ples with different degrees of disorder but with constant 
doping level. The systematic decrease of Tc with increas- 
ing disorder at fixed doping level indicates that out-of- 
plane disorder, rather than local modulation of the car- 
rier density, is a dominant factor that influences Tc . This 
is compatible with the recent STM-STS result by McEl- 
roy et alm^ on Bi2212 with excess oxygen atoms as a 
source of disorder. They conclude that disorder does not 
change the carrier density in the Cu02 plane locally. A 
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FIG. 4: (Color online) The variation of Tc as the function of 
(TA- The fitting line is obtained from the data of optimal doped 
Bi2212 and L-Bi2201 in our experiments. The plots of underdoped 
and overdoped Bi(Pb)2212 include the data from other reports 
(Refs^iaiaBiaB). The iusct shows the phase diagram of Bi-based 
HTSC with the axis of "disorder." 
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recent calculation shows that out-of-planc disorder, such 
as dopant oxygen atoms, could locally modulate the pair- 
ing interaction and work as impurities for ofF-diagonal 
scattering which gives rise to nanoscale modulation of 
A as well as to the suppression of the coherence peaksjiSi 
However, it is not clear whether strong disorder such that 
in Bi2201 may give arise to appreciable ZBC. 

In summary, a very broad gap distribution is observed 
by the STM and STS in the superconducting state of 
the disorder-controlled single-layer cuprate Bi2201. In- 
homogeneous gap structures were observed above and 
below Tc- As the out-of-plane disorder increases, the 
"pseudogap" region with larger gap and suppressed co- 
herence peaks increases while the superconducting region 
decreases. For the Bi-based cuprates including bilayer 
Bi2212, the degree of electronic inhomogencity measured 
by the width of gap distribution cta is found to show a lin- 
ear correlation with Tc when the doping level is kept near 
optimal. The result is suggestive of a possibility that the 



reduction of Tc with increasing cta is due to the reduction 
of pg. [In fact, it has been confirmed by the preliminary 
the muon relaxation rate (/iSR) measurement that ps in 
Gd-Bi2201 is lower than that in La-Bi2201(Refi)]. In 
this regard, the decrease of Tc with increasing cta is rem- 
iniscent of the relation between Tc and ps determined by 
/iSR, so called "Uemura-plot."— Thus, we speculate that 
the low Tc in Bi2201 is owing to strong disorder and Tc 
may be considerably improved if the disorder is reduced. 
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